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Spin-orbit interaction (SOI) in low-dimensional systems results in the fascinating property of 
spin-momentum locking. In a Rashba system the inversion symmetry normal to the plane of a 
two-dimensional (2D) electron gas is broken, generating a Fermi surface spin texture reminiscent of spin 
vortices of different radii which can be exploited in spin-based devices. Crucial for any application is the 
possibility to tune the momentum splitting through an external parameter. Here we show that in Pb 
quantum well states (QWS) the Rashba splitting depends on the Si substrate doping. Our results imply a 
doping dependence of the Schottky barrier which shifts the Si valence band relative to the QWS. A similar 
shift can be achieved by an external gate voltage or ultra-short laser pulses, opening up the possibility of 
terahertz spintronics. 

The Rashba effect 1 has been studied quantitatively on a variety of high-Z metal containing surfaces because of 
the accessibility of the surface states for spin- and angle-resolved photoemission (SARPES) experiments 2 . 
Although several surfaces show very large spin splittings, most of them have metallic substrates and are thus 
not suitable for the design of a spin-based devices 3 ' 4 due to the short-circuiting via the bulk states. Recently it was 
found that the surface of a topological insulator can also host Rashba-type spin split states that coexist with the 
non-trivial states 5 ' 6 . However, the origin of these states is likely due to adsorption processes and an increased 
interlayer spacing, and they are not susceptible to a small electric field 7 . One of the main reasons why it is difficult 
to control the Rashba effect in surface or interface states by an external electric field is that this field should act 
directly on the states themselves and thus fields in the order of 10 13 Vm" 1 are required for any substantial change 8 . 
For semiconductor heterostructures, on the other hand, the applied gate voltage indirectly affects the states via the 
symmetry of the confinement, and much lower voltages are required for a measurable change 9 " 11 . 

Here we propose that quantum well states (QWS) in an ultra-thin metal film grown on a semiconductor may 
offer an alternative way to realize a system where the Rashba effect can be manipulated by a small gate voltage. 
Measured momentum splittings are larger than in semiconductor heterostructures, indicating that much shorter 
channel lengths are required to achieve a similar spin-manipulation or to make a conversion between charge and 
spin currents. We focus on Pb films on Si(l 1 1) samples with different doping levels as a model system, because (i) 
Pb is a high-Z element (Z = 82) with large atomic SOI, (ii) smooth crystalline films can be grown on Si(l 1 1) with 
chemically sharp and well-defined interfaces 1213 , and (iii) the electric field in the depletion layer due to the ionized 
donors and acceptors can be varied by a gate voltage. 

Recently it was shown that QWS in thin Pb films on n-type Si(lll) show a Rashba-type spin splitting which 
builds up throughout the whole metal layer along the growth direction as a result of competing effects between the 
metal- substrate and metal- vacuum interfaces 14 . This is indicated by the absence of a strong dependence of the 
strength of the Rashba effect on the film thickness, and by the reversed spin helicity as compared to the surface 
state on Au(l 1 1) 15 . This is in contrast to other ultra-thin film systems where the spin splitting of QWS is induced 
by a high-Z metal substrate 1617 , or by hybridization gaps 18 . Later the sensitivity of the Rashba parameter to the 
film-substrate interface was determined in Pb QWS: replacing the Pb wetting layer by a Bi layer reduces the 
Rashba parameter by 60%, whereas QWS in Pb films grown on a Ag reconstructed Si substrate show no 
measurable Rashba-type spin splitting 19 . As will be explained below, these effects can be understood in terms 
of the local asymmetry of the QWS wave function with respect to Pb nuclei, which in turn is related to the 
particular confinement conditions. Here we will show that a similar effect can be produced by changing the donor 
concentration of the substrate which subsequently alters the Schottky barrier (SB) of the metal- semiconductor 
junction. 



SCIENTIFIC REPORTS | 3 : 1963 | DOI: 1 0.1 038/srep01 963 



1 



Results 

Figure 1(a) displays ARPES data from QWS formed in an N = 8 
monolayers (ML) thick Pb film deposited on lightly (upper panel) 
and heavily (lower panel) n-doped Si substrates [henceforth Pb/h- 
Si(lll) and Pb/1-Si(lll)]. The band dispersions are very similar to 
those obtained for Pb films grown on a moderately n-doped Si(lll) 
substrate [hereafter Pb/m-Si(lll)] 20 " 22 . Both band dispersions fea- 
ture an anomalously high effective mass (> 10 m e ), which is due to an 
increased in-plane lattice constant 23 . The films on both samples were 
measured with the COPHEE spectrometer 24 equipped with two Mott 
detectors, see Fig. 1(b), and were prepared by the same procedure. A 
change of the spin splitting is therefore directly related to the influ- 
ence of the donor concentration in the Si substrate, because the 
metal-vacuum boundaries are the same. Figure 2 shows SARPES data 
from a Pb QWS grown on the lightly n-doped Si(lll) measured for 
two energy distribution curves (EDCs) at (a) k|| = —0.11 A" 1 and (b) 
k|| = — 0.24 A~\ see also dashed arrows in the upper panel of Fig. 1 (a) 
which indicate the positions where the spin- resolved spectra were 
recorded. For the measurements at both momenta displayed in 
Figs. 2(c, d) the spin polarization data show typical up-down excur- 
sions in the tangential spin polarization components as expected for 
a Rashba system. 

A quantitative analysis of the spin-resolved data using a two-step 
fitting routine 25 reveals a larger spin splitting between the spin-up (f) 
and spin-down ( j) bands (AE : = e[ [k\\ ) — e\ (k\\ ) ) for the state with 
higher momentum. The corresponding fits are shown in Figs. 2 (a, c) 
for k|| = -0.11 A" 1 and in (b, d) for k|| = -0.24 A" 1 . We obtain a 
spin splitting of AE = -(22.2 ± 1.4) meV and -(33.8 ± 1.4) meV, 
respectively, with fully polarized bands. The increased spin splitting 
with increasing momentum is in full agreement with the Rashba 
model 1 where AE oc k||, and also seen in Figs. 2(g) and (h), which 
display raw data as spin-resolved EDCs projected onto the quantiza- 
tion axis at both momenta, calculated as 7^ = -7 tot (l +P tan ) with 
^tan = sign (Py) ^Jp 2 x + Pj. Figures 2(e) and (f) indicate the precise 

parallel momentum positions where the two EDCs were measured. 
For k\\ = -0.11 A" 1 the position was (k x , k y ) = (-0.086, 0.072) A" 1 
(Fig. 2(e)), for k\\ = -0.24 A" 1 it was (k x , k y ) = (-0.227, 0.072) A" 1 
(Fig. 2(f)). The strict Rashba-type spin-momentum locking requires 
the spins to be quantized in the direction perpendicular to the 
momenta, and the observed polarization along the x- direction is thus 
a consequence of the non-zero k y component. 



Figures 3(a-b) show SARPES data at k\ \ = -0.24 A" 1 for a Pb film 
of the same thickness prepared on the heavily n-doped Si(lll). The 
best self-consistent fit to the intensity and spin polarization data is 
here achieved with AE = — (48.0 ± 1.8) meV. A direct comparison of 
both systems at k\\ = —0.24 A" 1 reveals a larger energy splitting on 
the heavily n-doped substrate by almost 14 meV, which is well 
beyond the combined accuracy of 2 meV of the measurement and 
the fitting routine 25 . To illustrate that the larger AE of the QWS in Pb/ 
h-Si(lll) is not an artifact from our analysis procedure, we have 
performed a fit by keeping AE equal to the value found in Pb/1- 
Si(lll) (dashed lines in Fig. 3(a, b)). Although the total intensity 
fit is reasonably good, the fit to the polarization data obviously fails. 

From these and further data sets measured at different k\\ we can 
deduce the Rashba constants for both doping levels by plotting AE, 
obtained from the fitting, vs. k\ \ (Fig. 3(c)). Values for a RB obtained by 
using a RB = l/2-(dAE)/(dfc||) area RB , h = (0.11 ± 0.007) eVAforPb/ 
h-Si(lll), a RB , m = (0.076 ± 0.005) eVAfor Pb/m-Si(lll) 19 , anda RBa 
= (0.070 ± 0.011) eVA for Pb/1-Si(l 11). For Pb on a highly p- doped 
Si(lll) sample we obtain a RBa = (0.061 ± 0.006) eVA [Pb/p-Si(lll), 
N A = 1.18 • 10 16 cm" 3 , see SOM]. 

Discussion 

We now discuss the mechanism leading to the increase of the Rashba 
constant in Pb QWS on the n-doped Si(lll)with increasing donor 
concentration as shown in the right panel of Fig. 3(d). As demon- 
strated in Ref. 19 by the comparison of the Rashba effect in Pb QWS 
on a Pb and Bi terminated Si surface, a changed Rashba constant 
results from a modified charge density distribution in the Pb film 
which alters the local asymmetric features of the QWS wave function 
around the Pb cores. According to the phase accumulation model the 
charge density distribution of a QWS is controlled by the phase shifts 
at the metal-substrate and metal-vacuum interfaces 26 . Because the 
phase shift at the vacuum side is the same for all three investigated 
systems, we focus in the following on the metal-substrate interface 
with its phase shift given by ®s °c y/E — E 0 S (E — E 0 ) 13 where E is the 
energy of a QWS and E 0 is the valence band edge of Si at T and 0 is 
the Heaviside function. States with energies inside the gap of a sub- 
strate are truly confined, while states outside the gap can couple to 
substrate states leading to quantum well resonances with a consid- 
erable fraction of charge spilling into the substrate. In between these 
two cases the confinement, i.e. the degree of localization, changes as a 
function of E 0 . The goal of the following discussion is to show that the 




Figure 1 | (a) Band dispersion of a QWS in 8 ML thick Pb films on (upper panel) lightly n-doped and (lower panel) heavily n-doped Si( 111) substrates as 
measured by ARPES. (b) Schematic drawing of our SARPES experiment. 
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Figure 2 | Measured spin-integrated EDCs for Pb/1-Si( 1 1 1) (a) at ky = —0.1 1 A" 1 and at (b) k|| = —0.24 A" 1 with Voigt profiles obtained from the self- 
consistent two-step fitting routine, (c, d) Corresponding spin polarization data and fits of the x, y and z components. Insets show the in-plane polarization 
vectors obtained from the two-step fit. (e, f) Parts of constant energy surfaces illustrating the measurement positions of the EDCs in reciprocal space, (g, h) 
Spin-resolved intensity profiles projected on the quantization axis for the two measured momenta. 
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Figure 3 | Measured spin-integrated EDC for Pb/h-Si(lll) at k\\ = —0.24 A" 1 (a) with Voigt intensity profiles obtained from the two-step fit. 
(b) Corresponding polarization data and fits of the x, y and z components. The dashed lines are fits performed with AE = —33.8 meV. (c) Measured k- 
dependent energy splittings and linear fits (lines) to obtain ot RB for Pb/h-Si(lll) and Pb/1-Si( 1 1 1 ), respectively, (d) Measured oc RB vs. doping 
concentration (left panel) of the heavily p-doped Si(lll) and (right panel) of the differently n-doped Si(lll) substrates. 
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Figure 4 | (a,b) Energy diagrams of Pb/Si within the interface dipole model (not to scale) as a function of N D (upper panels), and capacitor-dipole 
model (lower panels), (c) Pb/Si in a capacitor to illustrate how £ ext is applied to the system (upper panel). Calculated Rashba constant of Pb/Si as a 
function of applied external E-field (lower panel), (d) Influence of £ ex t on the energetic position of the Si valence band edge, (e) Electric field in the parallel 
plate capacitor as a function of (reversed) gate voltage, (f) Charge density plots of QWS in Pb/Si with (solid line) and without (dashed line) E ext at 
k\\ = 0.25 A -1 , (g) Sketch of a spin-FET device operated with a gate voltage ( U G ) or by excitation with femto-second IR pulses (upper panel). Spin 
precession of % for electrons propagating along the x-direction in Pb films on l-Si(lll) and h-Si(lll) (lower panels). 



SCIENTIFIC REPORTS | 3 : 1963 | DOI: 1 0.1 038/srep01 963 



4 



energetic distance between the QWS and E 0 is sensitive to the 
Schottky barrier (0£ B ) and hence to the donor concentration because 
both quantities are related via: ® s n B +E 0 =E g where E g is the energy 
gap ofSi(lll). 

Our explanation is based on the interface dipole model for 
Schottky barrier formation, which takes the interaction of the metal 
and the substrate into account 27 . In this model the formation of polar 
Pb-Si bonds at the interface results in a charge transfer and the 
establishment of an interface dipole that directly contributes to 
the SB 27 : 

®? = ®M-Xs + eV M . (1) 

Here <I> M is the metal work function, Xs is me electron affinity of the 
semiconductor, and y int is the dipole induced voltage drop at the 
interface. We will argue that the size of V int is not only determined by 
the interface chemistry, but also by the doping concentration and the 
type of dopant (donor or acceptor). 

The upper panel of Fig. 4(a) shows the energy band diagram of Pb/ 
Si deduced from measuring the peak widths of QWSs and QWRs as 
suggested in Ref. 13, and from core-level shifts of the Si 2p spectrum 
in the interface region that revealed an accumulation of negative 
charge at Si atoms when the Pb-Si interface is formed. This is in-line 
with the predictions of the polarity of the Pb-Si bond using the 
Miedema electronegativity scale and with the study of Si 2p core- 
level shifts of Pb on n- and p-type Si that both revealed a shift toward 
lower binding energies 28 . According to Koopman's initial state model 
this means that the electric field generated by the interface dipole 
(Epb-si) points from Pb to Si independent of the type of dopant. The 
lower panel of Fig. 4(a) illustrates our model which explains the 
influence of the donor concentration on the band edge of Si(lll). 
It consists of the aforementioned interface dipole placed in a parallel 
plate capacitor (PC) that generates an electric field of strength 29 

|E PC | = yj 2eN D (U 0 + U G ) /ee 0 at the interface due to the space charge 
in the depletion layer and hence depends on the density of ionized 
donors (N D ). U 0 is the built-in potential across the Schottky barrier, 
U G is the gate voltage (for the moment U G = 0 V) and e= 11.9 the 
dielectric constant of Si. For Pb in contact with n-type (p-type) Si the 
polarity of the capacitor is such that the positive (negative) space 
charge is on the semiconductor side, while the negative (positive) 
charge of equal size is balanced on the Pb side. Hence the direction of 
the electric field of the dipole and of the capacitor are anti-parallel 
(parallel) to each other for n-type (p-type) Si. Now, increasing the 
donor concentration increases the electric field in the capacitor and 
lowers V int due to the screening of the interface dipole. Consequently, 
the SB decreases and the band edge of Si shifts to higher binding 
energies. 

Our finding of a reduction of the SB with increasing N D is con- 
sistent with previous results 30,31 . As already pointed out the increased 
energetic distance between the QWS and E 0 influences the confine- 
ment of the QWS via the metal-substrate phase shift such that the 
QWS is more strongly localized in the Pb film, i.e. less charge spills 
into the substrate. This affects the local charge densities of the QWS 
wave function with respect to the Pb cores which determine the spin 
splitting 19 . We now show how the charge density distribution in the 
Pb film changes as a function of the energetic position of the valence 
band edge using DFT calculations. These calculations were per- 
formed on a 10 ML Pb film on Si(lll) with the in-plane lattice 
constant of the Pb film commensurate with that of Si(lll) and with 
Pb/Si placed into a parallel plate capacitor producing external electric 
fields (E ext ) of various strengths with the positive bias at the Si side, 
see upper panel of Fig. 4(c). The DFT calculations reveal that increas- 
ing £ ext from 0 to 210 MVm" 1 shifts the Si valence band edge by 
~ 90 meV to higher binding energies, see Fig. 4(d), and simulta- 
neously the calculated Rashba constant increases by ~ 12%, see lower 
panel of Fig. 4(c). Notice that here E 0 shifts in the same direction as in 
the interface dipole model when the donor concentration is 



increased. As already speculated the increased Rashba constant is a 
result of an increased localization of the state in the Pb film, which is 
apparent from Fig. 4(f) where we show charge density distributions 
producing a a RB = 0.078 eVA at £ ext = 0 MVm" 1 and 0.087 eVA at 
210 MVm" 1 , respectively. This interpretation is also consistent with 
our results obtained from Pb QWS on heavily p-doped Si(l 1 1) where 
we find the smallest Rashba constant among the studied systems, see 
left panel of Fig. 3(d) and SOM. 

Having shown the sensitivity of the Rashba constant in Pb QWS to 
the donor concentration we now discuss the possibility to use Pb/Si 
as a candidate for future spintronic applications. Assuming ballistic 
spin transport through the Rashba channel we deduce the length (L) 
at which the injected spin - a coherent superposition of the two 
orthogonal Rashba spinors - precess by an angle of n while prop- 
agating along x, using L = n/(2k 0 ) where 2k 0 = 2m*oc RB /ti 2 is the 
characteristic Rashba-type momentum splitting deduced from the 
experiment. For electrons in Pb/1-Si(l 1 1) we obtain L/ = 1.71 nm 
fromkg,z = 0.092 A~\andinPb/h-Si(lll) = 1.08 nmfromko^ = 
0.146 A" 1 , see lower panels of Fig. 4(g). Such a difference in spin 
precession length can for example be used to design a spin-field effect 
transistor (spin-FET). For a spin-FET made of Pb/1-Si(l 1 1) with a 
fixed channel length of only 1.08 nm this means that changing the 
voltage drop at the interface by increasing the dipole screening via 
|Epc | gc y/Uo + U G it is possible to go from a low (A<p ~ n/2) to a high 
conducting state (A<p = n) by ramping U G from 0 to —12.8 V as 
simulated in Fig. 4(e). 

In a more general sense, the combination of the possibility to tune 
the Rashba effect via a gate voltage and the relatively large splitting 
make Pb QWS on Si an ideal system to study spin-transport effects, 
such as a varying spin accumulation due to the spin Hall effect 32 . In 
this respect it is important to note that an alternative way to vary the 
spin -splitting maybe achieved by optical means - thereby combining 
spin-based with electro- optical technologies. Recently it has been 
demonstrated that it is possible to change the SB and correspond- 
ingly E 0 in Pb/Si using ultra-short light pulses 33 . This would allow for 
a study of the time needed to build up or change the spin accumula- 
tion and indirectly the spin coherence time. Returning to our device 
example above, consider a spin-FET in the non-conducting state 
which can be optically driven. Illuminating the device with fs-pulses 
leads to a shift of E 0 such that k 0 changes and the conducting state is 
reached. For high-speed applications, the FET should respond 
quickly to variations of the trigger pulse. In Pb/Si a complete built- 
up of the E 0 shift is reached within 100 fs and equilibrium is recov- 
ered after 600 fs 33 - the non-conducting state. This yields a switching 
frequency in the terahertz (THz) regime. 

Our results demonstrate how the Rashba-type spin splitting of 
QWS in ultra-thin Pb films can be tuned effectively via the donor 
concentration of the Si(lll) substrate. We conjecture that the ener- 
getic position of the Si valence band edge is the decisive factor for this 
effect. Our results provide the first direct evidence of a doping 
dependent Schottky barrier. A similar shift of the VB maximum 
can also be achieved by other means, opening up the possibility to 
fine-tune the Rashba-type momentum splitting in metallic QWS via 
a gate voltage or by excitation with ultra-short laser pulses. 

Methods 

The ultra-thin Pb films were prepared in three steps consisting of (i) cleaning of the 
substrate by flashing above T > 1300 K to obtain the Si( 1 1 1)-(7 X 7) reconstruction, 

(ii) deposition of ~ 3 ML of Pb from a Knudsen cell and subsequent annealing at 
600 K for 60 sec. to form the dense (a/3 x y / 3)R30°(a) surface reconstruction 34 , and 

(iii) deposition of thin crystalline Pb films at T = 80 K. The films on both substrates 
were prepared by depositing the same amount of Pb with a rate of ~ 0.3 ML/min. The 
room temperature substrate resistance R and donor concentrations N D are for Pb/h- 
Si(lll): R = 4Qcm, N D = 1.14 -10 15 cm" 3 , for Pb/m-Si(lll): R = 53 Qcm, N D = 
8.38 -10 13 cm" 3 and for Pb/1-Si(lll): R = 95 Qcm, N D = 4.65 *10 13 cm" 3 . 

All measurements were performed with the COPHEE end- station at the Swiss 
Light Source of the Paul-Scherrer-Institute 24 at a base pressure below 3 X 10" 10 mbar. 
The sample temperature during the measurements was 80 K. The energy and angular 
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resolution of our spin-resolved (spin-integrated) ARPES measurements were set to 
100 (25) meV and ± 0.75° (0.5°) in respect of the low efficiency of Mott- scattering. In 
all measurements a photon energy of 24 eV was used because of the optimal 
photoemission cross-section for Pb QWS 35 . 

The ARPES spectrometer is equipped with two orthogonally mounted Mott 
detectors operated at an acceleration voltage of 40 kV 36 (see Fig. 1(b)). Each Mott 
detector measures asymmetries of the photoelectron current backscattered elastically 
off a polycrystalline gold foil, defined as A f = (I R — I L )/(I R + I L ), with / being the 
intensity measured in either right (R) or left (L) detector. The orthogonal arrange- 
ment of the two Mott detectors allows to measure asymmetries in A*, Ay and Ag 
directions of the Mott coordinate system (x,y,z). Three-dimensional spin polarization 
data is then obtained by dividing asymmetries by the Sherman function (S = 0.068), 
and by applying a rotation matrix (T) in order to transform the spin polarization 
components from the Mott reference frame into the sample coordinate system: 
Pi(E b ) = T(S- 1 A ] (£ & )) with i e (x,y,z). 

We employed density functional theory in the local density approximation 37 , using 
the full-potential linearized augmented plane-wave method in thin film geometry 38 
(for a program description see www.flapw.de). This allows us to expose the structure, 
consisting of a ten layer Pb film placed commensurately on six layers of Si(lll) 
substrate with a H -terminated backside, to an external electric field. This field is 
generated by effective condenser plates put in the vacuum at a distance of 5.3 A above 
and below the film structure. To simulate the Si-doping, we charge only the condenser 
plate below the Si and add the screening charge to the Si substrate. 
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